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Expansion of a wave function in a basis of eigenfunctions of a differential eigenvalue problem lies 
at the heart of the i?-matrix methods for both the Schrodinger and Dirac particles. A central issue 
that should be carefully analyzed when functional series are applied is their convergence. In the 
present paper, we study the properties of the eigenfunction expansions appearing in nonrelativistic 
and relativistic i?-matrix theories. In particular, we confirm the findings of Rosenthal [J. Phys. G: 
Nucl. Phys. 13, 491 (1987)] and Szmytkowski and Hinze [J. Phys. B: At. Mol. Opt. Phys. 29, 
761 (1996); J. Phys. A: Math. Gen. 29, 6125 (1996)] that in the most popular formulation of the 
i?-matrix theory for Dirac particles, the functional series fails to converge to a claimed limit. 
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I. INTRODUCTION 

Convergence of expansions of a one-component function into a series of eigenfunctions of a Sturm- 
Liouville problem was a subject of many studies. In some physical situations [l], 0], of particular interest 
are expansions of a function defined only on a finite and closed interval. The classical results on conver- 
gence of such series can be found, for instance, in [35] and for modern studies on this type of problems, 
including the equiconvergence method, the reader is referred to 0, E| and references therein. 

A similar problem for a two-component function was addressed by a number of mathematicians at 



rpc 

the beginning of the 20th century [9l-ll4| and later reviewed in numerous textbooks (see e.g. 
However, very few articles and textbooks deal with the development of an arbitrary function on a closed 
interval. Usually, either the expansion on an open interval is studied only or some additional conditions 
are imposed on the expanded function at the boundary points. To the best of the author's knowledge, 



the only classical paper discussing the general situation is the one by Birkhoff and Langer 12| ■ A more 



recent analysis of this kind of problems can be found, e.g. , in 15[ . A generalization of the equiconvergence 
Q 

method [7| to a vector case should be also possible. 

In the present paper, we apply the general results concerning convergence of eigenfunction expansions 
in the context of the i?-matrix theory of scattering processes. This theory was first developed for low- 



16- 18] for reviews on the 



energy collisions that could be described with the Schrodinger equation [l| (see 
subject). The i?-matrix theory for the Dirac equation [191 ] was formulated soon after the nonrelativistic 
one, with nuclear applications in view. Only later it was realized that the electron-atom collisions 
involving targets with large atomic numbers require a Dirac description, due to the increasing role of 



relativistic effects. The i?-matrix theory was reinvestigated in this context in 20j. 



The central idea in the formulation of the i?-matrix methods for scattering from spherically symmetric 



'Electronic address: julsta@ifae.es 



potentials is to divide the whole space into two regions, a finite reaction volume r < g and the outer 
region r > g, and to expand a wave function in the inner region in a series of eigenfunctions of the 
Hamiltonian governing the scattering process augmented, however, by artificial boundary conditions at 
the point r = g. This procedure allows one to express the i?-matrix as a limit r — > g~ of an infinite 
functional series. In both nonrelativistic and relativistic theories the critical issue, the convergence of the 
series on the boundary, was not properly analyzed by the originators and only presumed to hold. The 



convergence question was first recognized by Rosenthal 21], however his conclusions were incorrect. Later 



Szmytkowski and Hinze |22h24| ] realized that while the development of the solution in the Schrodinger 
formulation converges in the whole interval to an expanded function, the analogous series appearing in 
the relativistic case has a discontinuity at the crucial boundary point. In this way, in the most popular 
formulation of the method, the solution depends on the artificial boundary condition imposed on the 
basis functions. Taking this into account, Szmytkowski and Hinze developed the correct Dirac i?-matrix 



theory. Their conclusion caused much controversy [25| and was not widely recognized by the community 



261 ], We present here a theorem confirming their results as well as the general result on convergence 
obtained by Szmytkowski j^ . 

The paper is organized as follows. In Section [TTl we recall basic facts from both nonrelativistic and 
relativistic i?-matrix theories to highlight the problem of convergence appearing in both of them. In 
Section IIII1 we give the general convergence theorems concerning the eigenfunction expansions 0, fl"3 ) . 
The main result of the paper, a solution to the Dirac i?-matrix puzzle based on the theorem by Birkhoff 



and Langer [121 ]. can be found in Section HTT1 We finish the paper with conclusions and point out some 



open problems. 

II. EXPANSIONS APPEARING IN THE ii-MATRIX THEORIES 

The nonrelativistic and relativistic theories share many similarities, however in one essential point they 
are very different, i.e. the eigenfunction expansion of the solution of a nonrelativistic wave equation con- 
verges to a continuous function, whereas an analogous series in the relativistic theory has a discontinuity 



at the crucial boundary point 22j, l23l |27|. As a result, the relativistic i?-matrix is not appropriately 
expressed by a functional series. To highlight this difference, we shortly introduce both methods in 
a single-channel scattering from spherically symmetric potentials. The notation used in the following 
sections is based on the monograph on the i?-matrix methods in scattering 



A. Nonrelativistic ii-matrix theory 

A nonrelativistic elastic scattering process of spinless particles with mass m and energy E > from 
a spherically symmetric potential V(r) is governed by the stationary Schrodinger equation. We assume 
that the potential V(r) affects the particle only in a finite spherical volume of radius g centered at r = 0, 
denoted further by V e . Outside this (inner) region the particle is free and its wave function satisfies the 
free-Hamiltonian stationary Schrodinger equation. Obviously, the solution in the inner region must pass 
smoothly into the solution in the outer region. We denote by ^{E, r) the wave function being the solution 
of the respective Schrodinger equations in the inner and outer regions. 

Since the potential is spherically symmetric, it is enough to consider only the radial part of the function 
^(E,r) corresponding to the multiindex 7 = (I, mi). For a general function /(r), it is defined as 



F 7 (r) 



[ dfr 2 T 7 (r)/(r), ?=-, (1) 

Jilt r 



where T 7 (r) = (l/r)iT 7 (f), and Y 7 are normalized spherical harmonics defined as in [28[. We denote by 



P(E, r) a vector of radial functions of ^(E, r) with elements P 1 (E, r), and by D(E, r) a vector of radial 



2 



functions of f • V$(B,r) with elements D 7 (E,r). Let us assume that there exists a matrix R^(E,g) 
connecting P(E,r) and D(E,r) on the boundary of V g (which on the radial grid corresponds to r = g) 
in the following way: 

P(E, g) = R b (E, g)[D(E, g) - bP(E, g)}, (2) 

where b is an arbitrary square matrix. The above relation defines the R- matrix Rb(E, g). In what 
follows, it will be assumed that b is a diagonal, energy-independent, real matrix. Then the i?-matrix is 
also diagonal, and its elements will be denoted by (Rb) 77 = Rb-y- Finding the i?-matrix is equivalent to 



solving the scattering problem since Rb(E, g) is simply connected to the scattering matrix [18 1. 

To determine the eigenfunction expansion of the i?-matrix, we consider the radial part of the 
Schrodinger equation in the inner region V g : 

H 2 d 2 h 2 l(l + 1) , x \ , , r N 

We emphasize that here we do not make any restrictions on the function, except that it vanishes for r = 
as r l+1 . Our aim is to expand the unknown radial function P y (E,r) in the basis {P- (r)} generated by 
the same Hamiltonian, but augmented by the artificial boundary condition at r = g, that is 

h 2 d 2 , h 2 l(l + l) , . , _\ p(7 



2m dr 2 2mr 2 



+ V(r)-Ei P™(r)=0, r(=[0,g), (4) 



d 







= b 77 - 




r=g \ 


"5 



lim r-^ff J (r) = const, ^-Pp } (r) = ( b 77 + - ) P, (7) (g). (5) 



dr 

The set of eigenvalues Ei is countably infinite, and cigcnfunctions corresponding to different eigenvalues 
are orthogonal under the standard scalar product in C 2 ([0,g\). The formal expansion of an arbitrary 
function on [0, g] is 

P 1 {E,r)=Y, C ^ E ) P i l) ^ C ^ E )= drP^(r)P,(E,r), (6) 

where we assume that the functions {p/ 7 " 1 } are normalized to unity. The choice of the boundary conditions 
§5§ for the basis functions allows us to write the coefficient Ci (E) in a form that reveals the proportionality 
to D 1 (E, q) — b 77 P 7 (P, g) [compare to equation $2$}. In the case that we consider it holds that D~ f (E, r) = 
rd r {l/r)P 1 (E, r). Using equations ^ and (@| and the boundary conditions fulfilled by P/ 7 ^(r), we obtain 



1 

P,(E, r) = ^W^E, g) - b 77 P 7 (P, g)} £ -J_k|i>M( r ), r G [0, g). (7) 

i=0 1 



Taking the limit r — s- g on both sides leads to 

P 7 {E,g) = R b7 (E, g) [U 7 (E,g)- b 77 P 7 (E, g)} , (8) 

where 

Comparing equations ([SJ and © we see that equation Q defines a diagonal element of the P-matrix R b . 
One notices that the P-matrix is expressed as a continuous extension of a functional series to the point 
?' = g. The main question that we are going to answer in section 3 is: Can one interchange the symbols 
of limit and sum, and still obtain the same result? In other words, does the series on the right-hand side 
of equation ([9]) converge to a continuous function in [0, £>]? 



3 



B. Relativistic i?-matrix theory 



The relativistic description of an elastic scattering process for particles of spin | with rest mass m and 
total energy E {\E\ > mc 2 ) is governed by the stationary Dirac equation. Similarly as we have done 
in the nonrelativistic case, we assume that the potential V(r) vanishes outside the spherical volume V e 
bounded by a spherical shell, corresponding on the radial grid to r — g, and that the wave functions 
r) in the inner and outer regions pass smoothly one into the other on the boundary of V g . To define 
the relativistic i?-matrix, we fix the following notation: 







(10) 



where the multiindex 7 is defined as (k, rrij), with k £ Z\ {0} and rru = {— |k| + 1/2, — \k\ +3/2, . . . , |k| — 
1/2}, I = \k+ 1/2 1 — 1/2, and £l± Km Ar) are the spherical spinors [29j. We define two radial functions of 
a four-component vector f (r), denoted by the superscripts ±. For a fixed multiindex 7 they are given by 



F^(r)= [ dfr 2 r>Wt( r )f(r). 

J4-K 



(11) 



We denote by P(E, r) and Q(E, r) vectors of "+" and "— " radial functions of *f?(E, r), respectively, with 
elements P 7 (E,r) and Q~(E,r). Let us define the i?-matrix r£ (.E, g) connecting P(E,r) and Q(E,r) 
on the surface of V g in the following way: 



P(E,g) = R { b +) (E,g) 



2mc 



Q(E,g)-bP(E,g) 



(12) 



where b is some square matrix. Henceforward we will assume b to be diagonal, energy-independent and 
real. In this case the i?-matrix will be diagonal as well. 

To find the R- matrix, we consider the radial part of the Dirac equation in the internal region V g : 



mc + V(r) — E cH(-d/dr + n/r) \ / By(E,r) 
ch(d/dr + n/r) -mc 2 + V{r) - E I \ QJE, r) 



0, re [0,e). (13) 



Since the solution must fulfill some unknown boundary condition at the point r — g, we do not make 
any assumptions on the functions P 1 (E,r) and Q 7 (E,r), except for that they vanish for r — 0. We 
will expand these unknown functions in the basis generated by the eigenproblem consisting of the Dirac 
Hamiltonian of the previous equation and boundary conditions that, though unphysical, will allow us to 
develop the i?-matrix into a functional series. Let us consider the eigenproblem 



'mc 2 + V(r) - Ei cfi(-d/dr + n/r) \ ( P^\r) 
cH(d/dr + K/r) -mc 2 + V(r) - EJ Iq| 7) (r) 



0, re[0,f>], (14) 



lim r-»P^\r) = const, Qf\g) = {2mc/h)- 1 b 11 P^ ) (g), (15) 

1 7-0 



where v — I + 1 if V(0) — const, and v = k 2 — (aZ) 2 for the Coulomb potential (a is the fine-structure 
constant and Z the atomic number). The set of real eigenvalues {Ei} is infinitely countable, moreover, 
the eigenfunctions corresponding to different eigenvalues are orthogonal in the sense 

e dr(P^(r),Q^(r)) ( ^ ) =Af 2 6 tJ . (16) 

We will further assume the eigenfunctions to be normalized to unity; then the formal expansion of the 
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functions Py(E, r) and Q 7 (E, r) is given by 



Ry{E,r) 



V C(E)( P ^ {r) 



r G [0,g), 



(17) 



C t (E)= £dr(p^(r),Q^ 



(r) 



P y (E,r) 
Qj(E,r) 



(18) 



The coefficients Ci(E) can be written in the form revealing the connection to the i?-matrix. Using 
equations (fT5|) and (fH| and the boundary conditions fulfilled by P!f\r), we obtain 



Taking the limit r — >• 
where 



2m 



2mc 



Qj(E,g) - b 7J P 7 (E,g) 



Ei-E a (7) M 



re [0,£>). 



g on both sides we obtain for the upper component 



P^E, Q )=R$(E,q) 



2mc 



Q 7 (E,e) - b 77 F 7 (E,g) 



H 

ti 2 



P^(g)Prir) 



(19) 



(20) 



Comparing equation p^|) to equation p"2")) . we see that (f2"0")l defines the diagonal elements of the i?-matrix 



R y J>. Exactly as in the case of the nonrelativistic R- matrix (compare equation ©), the relativistic 



i?-matrix is expressed by a functional series whose convergence is directly related to the convergence 
properties of the series (|17[) . In the relativistic case the same question arises: Does the interchange of the 
limit and the infinite sum in equation still gives the same result? In other words, does the following 
identity hold: 



P^\g)P^\r) 7 
Ei-E 



E 



^ (7) (g)^ (7) (g) 

K-E 



(21) 



The expression on the right-hand side of the above equation is traditionally called the i?-matrix. However, 
in the next section we show that in the case of relativistic scattering it is not allowed to exchange the two 
operations. Therefore the "i?-matrix" as defined on the right-hand side of equation (|2~Tj) cannot connect 
the upper and lower component as in equation (|12p . 



III. CONVERGENCE THEOREMS 



In the previous section, we have reviewed the i?-matrix theories for the Schrodinger and Dirac particles. 
In both cases the i?-matrix has been defined as a limit of a certain eigenfunction expansion. In the 
nonrelativistic theory it is given by equation @, whereas in the relativistic theory by equation (I20[) . The 
problem of convergence of these two functional series is equivalent to convergence of expansions © and 
(TTTl) . respectively. Let us then recall the convergence theorems for these two problems starting with the 
one for the Schrodinger problem 

Theorem 1. (adapted from [S]) Consider the set of solutions of the following Sturm-Liouville problem: 



d2 



A, 



y„(x) = 0, 



a < x < b, 



y n (a) cos a + y' n (a) sin a = 
y n (b) cos f3 + y' n (b) sin f3 = 



(22) 
(23) 



5 



where q is assumed to be a real continuous function, and a, [3 £ [0, 271"). Suppose now that f is a real, con- 
tinuous function of bounded variation in the interval [a,b], then the expansion of f in the eigenf unctions 
of the eigenproblem \22\ ) + (2$) reads: 

oo 

f(x) = J2 c ny n (x), (24) 

with 

fb 



c n = / dxy n (x)f(x), (25) 

J a 

where the eigenf unctions y n are assumed to be normalized to unity in L 2 {[a, b]). The series f(x) converges 
uniformly to f(x) in the open interval (a,b). Moreover, for x — a and x — b the following relations are 
true: 

The above theorem can be applied to the nonrelativistic expansion, indicating that in equation © the 
limit can be exchanged with the infinite sum, giving 

R ( Eo) - fi2 y ^ 7) (^ (7) (g) (27] 
Rb 7 (S, g) - ^ ^ ■ (27) 

i=0 

In what follows, we present a theorem concerning the expansion of a two component function appearing 
in the Dirac i?-matrix definition. 

Theorem 2. (adapted from [12]) Consider the following boundary value problem: 

a<x<b, (28) 

(29) 

where A(x),B(x) are 2x2 matrices of functions continuous with their first derivatives, A(x) being 
diagonal, and W a ,Wb are constant 2x2 square matrices. Let us assume that: (i) the eigenvalues of 
A[x), denoted by *&j(x) (J — 1,2), are continuous functions fulfilling the following conditions for all 
x e [a, b] : 

Arg[i?i(.T) — *&2(x)] = const 

Argi9 i (a;) = const (30) 

Further, let us divide the complex plain of the parameter X into sectors in which the sign of the expressions 
5Rj = Re(A f dt'&jit)) (j = 1,2) is fixed. If one takes into account the conditions i30\) . there are two 
possibilities of dividing the complex plain of X. Both are visualized in figure\J\ where the aforementioned 
sectors are denoted by a^, k — 1, 2, . . .. For each sector we define the 2x2 matrices I and I with 
elements 

^— ky I 5ij, dtj <0 in sector at 
1 0, 3?j > in sector at 

(I% = { °' ^<0 in sector a fe . (31) 
I Sij , > in sector ak 



-[ 


u{x) 


Ax \ 


v v(x) 


w a 


' u(a) 




\ v{a) 
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ImA 



ImA 



a) 



\ 1 

0a(a:) \ 

\ <?2 


/ y iiti = o 


/ 

/ 

/ 

/ 

/ 

-1 


1 

5R 2 = ^ 

(J 4 



ReA 



b) 



1 


/ / 


/ 02 


— ~^s/ / 




i 






SRj = / / 




/ / 




/ / 








/ 





ReA 



FIG. 1: Two possible ways of dividing a complex plain of the parameter A into sectors in which the value 9?j has 
a fixed sign. In figure a) the two eigenvalues fulfill the condition Arg[#i(a;)] — Arg[$2(:r)] 7^ tv while in b) the 
two eigenvalues have phases differing by ir. The continuous halflines mark the rays on which the values of r t}j{x) 
lie, whereas the dashed lines — the rays on which the functions 3lj are equal to zero. The dashed lines are the 
borders of the sectors. 



(ii) Let the matrices W a and Wb be such that in each sector a k the following matrix is invertible: 



n k = w a i k + w b i k . 



(32) 



Then: 1. the set of eigenvalues A„ and respective normalized eigenf unctions (u n ,v n ) T of the problem 
) + i29\) is infinitely countable; the orthonormality relation is the following: 



dx{u n {x),v n {x))A{x) 

\ v m {x) 



(33) 



where (u n (x), v n (xf\ is the eigenvector of the adjoint boundary value problem )3a l. corresponding to the 
eigenvalue A„. 2. the development of any two- component function F = (Fi,F2) T , real and continuous 
with the first derivative in the interval [a, b], in the eigenf unctions of the boundary problem H28\) + Wtyl is 
given by 



F(x)= C " 



n— — 00 



u n (x) \ 
v„{x) I ' 



with 



L 11 



dx(u n (x),v n (x))A(x) ( *^ ' X . 



The expansion \3J$ has the following propertie 



F(x) = F(x), for a < x <b, 
F(a) = H a F(a) + J a F(b), 
F(b) = H b F(a) + J b F(b), 



(34) 



(35) 



(36a) 
(36b) 
(36c) 



where the 2x2 matrices H a , J Q , Hb, Jb are fully determined by the matrix A(x) and boundary conditions, 



7 



and given by the expressions 



Ja 
H b 

Jb 



ii 1+ 5:[-£TV"'j' 

k 

Er wfe- 
L 2tt 

k 

£[- 

it 



27T fc 



(37a) 
(37b) 
(37c) 
(37d) 



In the above, the parameter uik is an angle between the boundary rays of a sector ( dashed lines in Fig. 
Q]) and J\- is the invertible matrix defined in equation 



The above theorem is the special case of a more general eigenfunction expansion problem considered 
in 12j. The reader may find there the proof of the above facts. We will present here in more details 



the situation directly applicable to the i?-matrix expansion. Let us than state and prove the following 
corollary. 



Corollary 1. Consider the boundary value problem 

I p(x) — Xp(x) — d/dx + t(x) 
\ d/dx + t(x) q(x) — \p(x) 



cos a sin a 




/(«) 
9(a) 



f(x) 

o o" 

cos p sin /3 



= 0, a < x < 6, 



f(b) 
9(b) 



0, 



(38) 
(39) 



with p, q, t, p being real junctions continuous with first derivatives, p(x) > for all x € [a, 6] and a, /3 € 
[0,27r). Then the set of eigenvalues A„ and eigen] 'unctions (f n ,9n) T is discrete. The expansion of a 
two- component function F — (F%, F%) T , continuous with first derivative in [a,b], in the set {(f n , 9n) T } is 
given by 



oo 

F(x)= C " 
n=— oo 
rb 



fn(x) 
9n(x) 



C n = dx p(x){f n (x),g n (x)) 



Fi(x) 
Fi(x) 



and has the following properties: 

F(x) = 

F(a) = 

F(b) = 



F(x), for a < x < b, 

1 / 1— cos 2a —sin 2a 

2 I — sin 2a 1 + cos 2a 

1 / 1 -cos2^ -sin2/3 

2 I - sin 2(3 1 + cos 2/3 



F(a), 
F(b). 



(40) 

(41a) 
(41b) 

(41c) 



Proof. To prove the corollary, let us rewrite equation (f3"5|) in such form that the results from [12| , recalled 
in Theorem^ apply directly, i.e. we would like to have the differential equation and boundary conditions 
in the form (|28[) . To achieve this, we multiply equation (|38[) on the left-hand side by the unitary matrices 
U and U given by 



U 



- r 1 

V2 \ 1 i 



U 



i 
-i 



(42) 



<S 



obtaining 




ip(x) 
-ip(x) 



B(x) 




(43) 




where u = (if + g)/V^ and v = (/ + ig)/y/2, the diagonal matrix on the right-hand side corresponds to 
the matrix A, and the matrix _B contains the functions p,q,t, and fulfills the assumptions of Theorem [5] 
At the same time, we have to adjust the boundary conditions to the functions u and v. These become 

= 0. (44) 

Comparing equation (|44[) to (|29[) . we can see that the 2x2 matrices on the left-hand side can be identified 
with W a and Wb, respectively. Let us then check if equation (j4"3"]l with boundary conditions (|4"4"]) satisfy 
the assumptions (i) and (ii) of Theorem [2] 

First, we will find the division of the complex plain of the parameter A into sectors, as described in 
Theorem^ The matrix A(x) has two complex eigenvalues i)i(x) = ip(x) and ^(x) = ^p{ x )- Note that 
since p(x) is strictly greater than zero and continuous for all x € [a, b], they satisfy the condition (i) of 
Theorem [5] Note, in particular, that Arg[i?i(x)] — Arg[i?2(x)] = tt, which leads to the division of the 
complex plain of the type shown in Figured}), i.e. we have the following two sectors: 

o-i = {A : sfti < OA5R 2 > 0} = {A : ImA > 0} (45) 

and 

er 2 = {A : 3?i > A 9f? 2 < 0} = {A : ImA < 0}. (46) 
This implies that the matrices defined in (I3T1) are 

I 2 = 

Now the condition ([3^)1 can be checked easily. Let us write out explicitly the matrices Sli and for W a 
and Wb as defined in (I44j) . since we are going to use them in further calculations: 

/ gia e -i/3 \ / _j i/3 j -ia \ 

! - ( r "^[^ ■-*)• (48) 

Clearly, they are invertiblc for all values of a,/3 € H, so the condition (ii) from Theoremr2]is fulfilled. 

We have checked that all assumptions of Theorem [5] are satisfied, therefore, the set {{u n ,v n ) T } of 
eigenfunctions of the problem (l43|) + (j44"]) [and at the same time the set {(f n , 9n) T } of eigenfunctions of the 
problem (|3"5|) + (|59")) ] is countably infinite. Moreover, part 2 of the theorem applies to the expansion (|40l) . 
To show explicitly that formulas (|41a[) - (|41c[) are valid, we premultiply equation (j4"0)) by y/2U obtaining the 
series of (u n , v n ) T representing a complex two-component function G = (Gi, G*2) T = (F 2 +iFi, F\ +iF 2 ) T : 

G(x)= V C n ( Un{x A. (49) 






Although part 2 of the theorem was formulated for real functions, the generalization to complex functions 
is straightforward if the real and imaginary part are considered separately. Therefore, we will proceed 
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with the complex function G. Notice that coefficients C n can be written as 



C n = f &xp(x){f n {x),g n (x))uWUU 



dx(u n (x),v n (x))A(x) > ( 5 °) 

where (u n (x), v n (x)) T = (l/2)U(u n ,v n ) T is the solution of the eigenproblem adjoint to (pf*3")) + (|4"4"]) . The 
obtained formula for coefficients is in agreement with equation (j3"5)l . Consequently, the series ([^Ujl modified 



with U converges to expressions (|36a[) . (|36bl) . and (|36cl) . Let us then determine the matrices H a , J a ,Hf,, Jb 
in this particular case. The angles co\ and w 2 are equal to n, so exploiting equations (|37al) - (|37d[) we obtain 

Ha = l \-X l \ )' Hb= (oo)' 

_(0 \ _J 1 -ie-^\ (r>1) 

V J ' 2 JUf 3 1 ) ' 

Taking into account the above results and equations (|36b[) and (|36c[) , we obtain the following expressions 
for the sum of series ([4"9")) at the points x = a and x = b: 

G(a) = H a G(a), (5(6) = J 6 G(6). (52) 

We recover the formulas for the original function F premultiplying equations (|52l) with (l/v2)J7^: 

F(a) = U^H a UF(a), F(b) = UU b UF{b). (53) 

Finally, we insert the matrices (|S"Tj) into the above equations and obtain slightly modified formulas (I41b[) 
and BTc|) : 



1 / 2 sin 2 a -sin2a\ / F x (a) 

2 y-sin2a 2 cos 2 a / \F 2 (a) 

1 / 2 sin 2 /? -sin2/3 \ I Fi(b) 

-sin 2,3 2 cos 2 /3 J 1 F 2 (6) 



— , . _L I j£j Dill Lt Dill ^Lt I I 1 1 1 HI | . . 



Summarizing, in the open interval (a, b) the series (|40|) converges to the function F(x), whereas at the 
boundary points x = a and x = b, the sum of the expansion is given by (I54a[) and (|54b[) . respectively. 
This finishes the proof of Corollary [TJ □ 

The Corollary reveals that the series (j4"0)) converges to the function F(x) in the whole interval [a, b] if 
and only if the following two equalities hold simultaneously: 



1 / 2sin 2 a -sin2a\ / F 1 (a)\ _ ( F x {a) 

2 ^-sin2a 2cos 2 a J \F 2 {a) J ~ \ f 2 {o) 

1 / 2 sin 2 ,3 -sin 2/3 \ / Fi(6)\ _ / F x (b) 

2 1 - sin 2/3 2 cos 2 /3 J 1 F 2 (&) J ~ \F 2 (b) 



(55) 
(56) 



which is equivalent to the condition 

F x (a)cosa + F 2 (a)sma = 0, F x (b) cos/3 + F 2 (b) sin/3 = 0. (57) 

One recognizes in this formulas the boundary conditions (1391) . Consequently, the sum of the series is 
continuous if and only if a function to be expanded fulfills the same boundary conditions as the basis 
functions. 
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We demonstrated in this section that the properties of the developments into eigensolutions of first- 
order differential systems (Theorem 0) and the expansions in the eigenfunctions of the Sturm-Liouville 
problem are dramatically different. This fact, not realized by the originators of the relativistic i?-matrix 
method, has far-reaching consequences for the theory. Although the assumptions in Corollary [1] are 
restrictive, i.e. the functions p, q, t, p, being elements of the matrices appearing in the eigenproblem, are 
assumed to be continuous with their first derivative, its conclusion still can be applied to expansion (|17|) 
at the point r = g. This can be done, because the boundary conditions are separated and, consequently, 
the convergence of the eigenfunction expansion at the point r = g is independent of the behaviour at the 
point r = 0. Taking this into account, we immediately see that the expansion (fl~7| for r = g does not, 
in general, converge to the solution of (fT3)l . It does only if the functions P 7 (E,r) and Q 1 (E,r) satisfy 
the second of the boundary conditions (1151) . This, however, cannot be assumed since this particular 
condition does not have any physical meaning and is chosen in this way only to obtain the expansion 
of the i?-matrix. As a result, the commonly used definition of the i?-matrix contains an error, since in 
general it holds that 

- g^V) - 

i— — oo i— — oo 

The way to correct this mistake was found by Szmytkowski and Hinze for a general multichannel case 



18l |22h24| . They introduced the correction which should by subtracted from the common and faulty 



expression for the i?-matrix in order to obtain the correct one. 



IV. CONCLUSION 

Summarizing, we have provided theorems concerning the convergence of eigenfunction expansions of a 
two-component function into eigenfunctions of a Dirac operator on a finite closed interval augmented by 
separated boundary conditions. In particular, we have shown that such expansions have discontinuities at 
the boundary if the expanded function does not fulfill the same boundary conditions as the basis functions. 



This confirms the result of Szmytkowski 27] and has far-reaching consequences for the relativistic R- 
matrix method. Moreover, the fact that the functional series does not, in general, converge to a continuous 
function in the closed interval may affect the rate of convergence and cause the Gibbs-like phenomenon 
3C, 31| to occur. 



The issue left as an open problem is the proof of convergence of (]40fl in the case when the functions 
p(x),q(x),t(x) have a singularity at one of the boundary points, e.g. at x = a. This is directly related 
to convergence of expansion (|17l) since the functions have discontinuities for r = 0. However, due to the 
separated character of the boundary conditions (|39f and the fact that both the expanded function and 
the basis functions vanish at r = 0, the conclusions (I41ap and (|41c|) of Corollary [1] should hold even in 
this case. 
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